Common Fermi-liquid origin of T 2 resistivity and superconductivity in n-type SrTiOs 
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A detailed analysis is given of the T 2 term in the resistivity observed in electron-doped SrTi03. 
Novel bandstructure data are presented, which provide values for the bare mass, density of states, 
and plasma frequency of the quasiparticles as a function of doping. It is shown that these values 
are renormalized by approximately a factor 2 due to electron-phonon interaction. It is argued that 
the quasiparticles are in the anti-adiabatic limit with respect to electron-phonon interaction. The 
condition of anti-adiabatic coupling renders the interaction mediated through phonons effectively 
non-retarded. We apply Fermi-liquid theory developed in the 70's for the T 2 term in the resistivity 
of common metals, and combine this with expressions for T c and with the Brinkman-Platzman-Rice 
(BPR) sum-rule to obtain Landau parameters of n-type SrTiOs. These parameters are compara- 
ble to those of liquid 3 He, indicating interesting parallels between these Fermi-liquids despite the 
differences between the composite fermions from which they are formed. 
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I. INTRODUCTION 

SrTiOa is a semiconductor which, when doped with a 
low density of electrons, becomes a good conductor with 
relatively high mobility and strong temperature depen- 
dence of the electrical resistivity and the infrared optical 
conductivity. At low temperatures the material becomes 
superconducting^ with a maximum reported T c of 1.2 
although superconductivity is usually reported below 0.7 
K with a dome-shaped doping dependence of Tj^l Su- 
perconductivity is also observed below 0.3 K in the two- 
dimensional electron-gas formed at the interface between 
SrTi0 3 and LaAlC# where the carrier-concentration de- 
pendence of T c has also a dome shaped The DC resis- 
tivity below 100 K has a T 2 temperature dependence, 
which has been attributed to electron-electron scatter- 
ing by some groups^^. On the other hand, resistivity 
of the form p(T) oc 1/sinh 2 (ujq/2T), which is almost 
T 2 -like, was found in Lai_ K Ca x MnOaf^ and in doped 
LaTiOjp' in accordance with the expected behavior of 
small polarons^. However, n-type SrTiC>3 appears to 
be described well by the model of large polarons with a 
Froehlig-type electron-phonon interaction^. Conditions 
in this material are therefore rather remote from those 
addressed by the small-polaron modeP^, and the ques- 
tion as to why the T 2 behaviour dominates up to high 
temperature remains as yet open. 

The resistivity near absolute zero has been known 
to be of the form p = AT 2 in platinurrfS] anc j other 
transition metal elements^^, with A ranging from 
2.5 • 10- 6 /i^cmK~ 2 (osmium) to A = 3 ■ 10" 5 /x^cmK- 2 
(palladium). M. J. Rice has expl ained these observations 
in terms of the Baber mechanis m 1 20 1 2 1 1 T 2 resistivity was 
subsequently observed in the alkali metals (s ee Re f.l2"2"lfor 
a review), with A = 3 • lO-V^cmK" 2 for Lp^H and an 
order of magnitude smaller values for K and N a 22 | 25 | 2 Sl . 
Based on the assumption that the Coulomb repulsion 



is the only interaction between electrons, Lawrence and 
Wilkins?3 calculated values in the range from 10 -8 to 
10~ 10 /iS!cmK~ 2 for the alkali-metals. MacDonald ob- 
tained similar values, and showed that the dominant con- 
tribution to the T 2 term in the re sistivity results from 
phonon-mediated interaction ;] 28 ! 29 !. A value several or- 
ders of magnitude higher, A = 0.02/i£lcmK~ 2 , was ob- 
served for stoichiometric TiSjp -', and the resistivity of 
Ti 1+2 .S2 as a function of carrier concentration was ob- 
served to follow the relation n~ 5 / 3 T 2 in agreement with 
the theoretical expressions in Ref. 27. 

In 1968 M. J. Rice pointed outP, that the coeffi- 
cient A should vary predominantly as the square of the 
linear electronic specific heat coefficient 7; in particu- 
lar he showed that the experimental data of elemen- 
tal 3d, Ad and 5d transition metals satisfy the relation 
A/j 2 — A • 10~ 7 pil cm (mole K /mJ) 2 . Heavy fermion 
compounds are characterized by very large values of A 
and 7. Kadowaki and Woods 31 summarized the situation 
by showing that A/7 2 in this group of materials is a fac- 
tor ~ 25 larger than in aforementioned data of elemental 
transition metals. Accordin g to the theory of electron- 
electron scatteringSSEIEZEIl th e ratio A/7 2 contains in- 
deed several non-universal factors, including the square of 
the strength of the effective electron-electron interaction. 
Since in general the interactions differ in nature from one 
group of materials to another, the same values of A/7 2 
are only expected within a particular group. The carrier 
density constitutes another non-universal factor, which is 
particularly significant for doped semi-conductors in view 
of their tunable carrier density. Hussey 9 -' proposed there- 
fore a re-scaling of the Kadowaki- Woods plot to account 
for, among other factors, variations in carrier density, and 
demonstrated that this notion is supported by the strong 
doping dependence of A in hole-doped LaTi03. 

Here we return to the possibility that the T 2 resistiv- 
ity in n-type SrTiOs could be a consequence of a quasi- 
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non-retarded interaction between dressed quasiparticles. 
The A-coefficents of SrTii-^Nb^Os, a few examples of 
which are listed in table [TTJ arc large. Since, as has been 
demonstrated by Thompson 30 , A oc rt~ 5 / 3 , this is a nat- 
ural consequence of the low carrier density. For exam- 
ple, SrTio.98Nbo.02O3 has a carrier density n — 3.4 • 10 20 
cm -1 , while lithium n = 4.7 • 10 22 cm -1 . If we as- 
sume that everything else is the same for these two ma- 
terials, the A coefficient of SrTio.9sNbo.02O3 should be 
4000 times larger than the one of Li. In reality they 
differ by a factor 8000, hence from this perspective the 
strength of the quasiparticle-quasiparticle scattering in 
SrTii-^Nb^Os is not drastically different from that in 
lithium. 

An obvious source of interaction in doped SrTi03 is 
provided by the overlap of the screening clouds surround- 
ing the electrons provided by the interaction with the lat- 
tice. The main phonons involved in this screening are op- 
tical ones, with the important consequence that their en- 
ergy exceeds the Fermi energy for the doping levels where 
superconductivity is observed. The polaron-polaron in- 
teraction mediated by these phonons is then effectively 
non-retarded, an unconventional aspect which we con- 
sider to be crucial for the observed T 2 dependence of 
the relaxation rate. The effective electron-electron inter- 
actions can also lead to the formation of Cooper pairs. 
Based on our analysis of the T 2 relaxation rate and of 
the superconducting transition temperatures we obtain 
an interaction of weak to moderate strength, making im- 
plausible scenarios where a substantial fraction of the 
charge carriers is paired in the normal state. 



II. TRANSPORT PROPERTIES 

In Fig.[T]the transport data of SrTi 1 _ 2 ,Nb 2; 03 with dif- 
ferent carrier concentrations are shown as a function of 
temperature^. Hall data are presented as Rh,o/ Rh(T), 
where Rh,o represents the zero temperature limit, for 
which the Hall charge carrier densities per unit cell xr — 
-a?/(eR H ,o) are 0.105%, 0.196%, 0.870% and 2.00%, 
which is within 12 % of aforementioned Nb concentra- 
tions specified by the supplier. At 4 K we observe fairly 
high mobilities in the range from 400 to 6000 cm 2 /Vs, 
which drop gradually as a function of increasing temper- 
ature to approximately 6 cm 2 /Vs at room temperature. 
These high mobilities at cryogenic temperatures are the 
first indication that n-type SrTiC>3 is a clean Fermi-liquid 
of mobile charge carriers. Concentrating now on the tem- 
perature dependent properties, we take a closer look at 
the inverse Hall constants. First of all we notice that 
the sensitivity to temperature changes diminishes for in- 
creasing carrier concentrations. We consider the possi- 
bility that the system has multiple electron-type bands. 
The effective Hall density n# = —e/Rff of a multiband- 
band system with carrier density n and fractional occu- 
pation of the j'th band Xj with mobily fij is given by 
the expression nn/n — x j^j) 2 x jf 1 'j) — 1- The 
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FIG. 1: Temperature dependence of the resistivity, the inverse 
Hall constant and the mobility of SrTii-^NbzOa for different 
carrier concentrations. 



limiting case nnjn = 1 occurs when only one band is 
occupied, or/and if fij is independent of j. In all other 
cases nn/n < 1. The temperature dependence is well 
illustrated by the case where at T = only one band is 
occupied (afi = 1). Increasing temperature makes Xj > 
for j > 2 and xi < 1, consequently nn/n is reduced. 
When two or more bands are already occupied at T = 0, 
the relative change in occupation number as a function of 
temperature is weaker and consequently nn will be less 
temperature dependent. In Ref. [34] a weak temperature 
dependent decrease of the Drude spectral weight, w 2 , was 
reported for temperatures higher than 100 K. A grad- 
ual temperature induced transfer of part of the electrons 
to states with a higher effective mass (and consequently 
lower mobility) then provides a natural explanation for 
both phenomena: The temperature induces a decrease of 
uip because it is inversely proportional to the mass, and 
an increase of Rh ■ 

The resistivities have a small residual component. The 
values for po were determined by fitting the data below 
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15 K to a constant plus a power law, and these values of 
Po are used in the remainder of the analysis. The inset 
of Fig. [T] shows the resistivity, from which the residual 
component has been subtracted, on a double log scale, 
indicating a power law like increase as a function of tem- 
perature. For further analysis it is useful to convert the 
resistivities to relaxation rates using the expression 
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For ujp we substitute the values measured with time- 
domain infrared spectroscopy on the same set of 
samples 34 . The residual relaxation rate turns out to be 
proportional to x. Since x is just the density of Nb 4+ 
ions, and these ions act as scattering centers, this (near) 
proportonality of scatttering to Nb-concentration is rea- 
sonable. In the following discussion we will focus on the 
behavior of h/r below 100 K, where both and Rh are 
independent of temperature. The results of least-square 
fitting the relation h/r — H/tq = a n (T/100K) v , summa- 
rized in table[lj clearly demonstrate that the temperature 
dependence of the resistivity up to 100 K follows closely 
a T 2 power law. With this in mind we fitted a 2 in the 
expression h/r — h/ro — a 2 T 2 , which values are listed in 
table [Tl] and the corresponding fits are displayed together 
with h/r in Fig. [2j Attempts to improve the fit in the 4- 
100 K range by adding a T 3 term decreases x 2 , and affects 
a 2 somewhat. However, the pref actor of the T 3 term is 
negative for x=0.009, hence a T 3 term below 100 K gives 
unphysical results and should be dropped, with the only 
possible exception the x=0.02 sample. The difference be- 
tween the data and the fit is constant upto 100 K, and 
grows rapidly at higher temperature, indicating that an 
additional component to the resistivity becomes active 
at that temperature. Such behavior is consistent with 
aforementioned interpretation of Hall data and spectral 
weight data, namely if electrons are transferred to lower 
mobility states, the resistivity will deflect upward from 
the trend observed at lower temperatures. We will return 
to this issue in the discussion of the mean free path in 
section Ivl 
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FIG. 2: Temperature dependent relaxation rates of 
SrTii-^Nb^Oa for 4 different carrier concentrations, using 
the relation p(T) = Atvuj~ 2 t~ 1 . Plasma-frequencies, lj p , are 
obtained from the Drude spectral weight measured with in- 
frared spectroscopy^ and listed in table [il] Inset: Differ- 
ence between experimental data and fitted curve, x(T) = 
h/r(T) — h/rfit(T), demonstrating upward departure from 
T 2 behavior of the resistivity above 100 K. 
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TABLE II: First column: Hall numbers rounded off to one 
significant digit, used in Figs. [2] and [7] to label the samples. 
Second and third columns: Fitting parameters and variance 
of the temperature dependent relaxation rate to a T 2 law. 
Fitting curves corresponding to 02 are compared to the ex- 
perimental data in Fig. [2] Columns 4 to 6: Fitting parameters 
and variance of the temperature dependent relaxation rate to 
a T 2 + T s dependence. Column 7: The A coefficient in the 
relation p = po + AT 2 . The values of h/ro are those of tablejl] 
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TABLE I: First column: Nominal doping. Second column: 
Hall number in the zero temperature limit. Third column: 
Residual resistivity. Fourth column: Drude plasma frequency. 
Column 5: Residual relaxation rate. Columns 6 to 8: Fitting 
parameters of the temperature dependent relaxation rate fit- 
ted to a power law and corresponding variance. 



III. BAND STRUCTURE 

SrTiC>3 has a cubic crystal structure at room tem- 
perature, which becomes tetragonal below a structural 
phase transition at 105 K. A 3 eV gap separate s the 
filled oxygen 2p bands from the empty Ti 3d bands^iM'. 
In Refs. 13 34 37 we compared experiments to novel ab 
initio band calculations, the details of which have not 
been presented in the literature. In the present article 
we make again extensive use of the same new ab initio 
data. Since there are some differences compared to pre- 
viously published bandstructure calculations, the new ab 
initio calculations are presented here in some detail. 

First principles calculations were performed using the 
Linear Augmented Plane Wave method as implemented 
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in the WIEN2k code^ and the generalized gradient ap- 
proximation for the exchange-correlation potential in the 
form proposed by Perdew and coworkers^! (See Ap- 
pendix [A]) A detailed view of the bandstructure around 
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FIG. 3: Band dispersion of the lowest unoccupied bands of 
SrTi03 in the low temperature tetragonal phase. The direc- 
tions in momentum space are labeled according to the high 
temperature cubic Brillouin zone, so that [1, 0, 0] corresponds 
to momentum along the Ti-O bond direction. The rightmost 
panel indicates the position of the Fermi energy as a function 
of carrier concentration. x cl — 4.0 ■ 10 -5 and x C 2 = 2.6 ■ 10~ 3 
are critical carrier concentrations where the Fermi energy en- 
ters the second and the third band. 




FIG. 4: Fermi surface of the high-temperature cubic phase at 
2% doping, showing the large anisotropy of the lowest band. 
At the critical doping x c = 0.097 a topological transition takes 
place where the Fermi surfaces open up along the three axis. 

the zone center is shown in Fig. [3] In this limited region of 
k— space the band structure can be effectively described 



by a tight-binding model within the ti g manifold of the 
Ti-3d states. The main aspects of the band structure 
are described by Bloch- waves of d xy , d yz or d zx charac- 
ter, each of which has two directions of strong dispersion 
(k x and k y for the d xy orbital etc.) and one slowly dis- 
persing direction orthogonal to these. The result is a set 
of three degenerate bands. The Fermi surface consists 
of three interpenetrating ellipsoids centered at the zone 
center, with the ellipsoids oriented along the x,y and z 
axis of the reciprocal lattice of the cubic crystal structure. 
This zone-center degeneracy is however lifted by the spin- 
orbit interaction. In the cubic phase this results in two 
degenerate spin-orbit doublets at the lowest energy, and 
an additional doublet at 29.2 meV higher energy. This 
splitting equals 1.5£, where £ = 18.8 meV is the spin- 
orbit parameter, somewhat smaller than £ = 25 meV 
used by Mattheiss. In the low temperature tetragonal 
phase the crystal field D — 2.2 meV lifts the degener- 
acy between the two doublets causing a splitting of 4.3 
meV. The result is the following set of bands having their 
minimum energy at the zone center: The lowest "heavy 
electron" band consists of states carrying angular mo- 
mentum rrij = ±3/2(1 — 6), where S cx D 2 /£ 2 . While the 
band disperses upward rather sharply at the zone center, 
it is deflected downward at \k\ rs 0.1 /a for momentum 
along the Ti-0 bond. The second band is a "light elec- 
tron" band, which becomes occupied at the critical car- 
rier concentration x c \ = 4.0 ■ 10~ 5 . Its dispersion is to 
a good approximation an isotropic parabola, and these 
bands have the peculiarity that the gyromagnetic factor 
gj = due a compensation of orbital (<?; = 1,to; = ±1) 
and spin magnetic moment (g s — 2,m s = ±1/2). The 
third band is also a light electron band which becomes oc- 
cupied at the critical carrier concentration x C 2 = 2.6-10 -3 
(n=4.4 ■ 10 19 cm~ 3 ). An experimental indication for this 
critical carrier concentration comes form the observation 
by Binnig et alW of an additional superconducting gap of 
smaller size than the main gap for doping concentrations 
in excess of 5 • 10 19 cm -3 , using tunneling spectroscopy. 

The most signicifant differences between the results 
presented here and Matheiss' result^! are the much 
smaller crystal field parameter D = 2.2 meV obtained 
here as compared to D = — 33 meV obtained from a tight- 
binding fit to Matheiss' bands, and the fact that the sign 
is opposite. The resulting Fermi surface of the lowest 
band is therefore quite different: In the present calcu- 
lation it is in fact similar to Fermi-surface of the cubic 
phase shown in Fig. E] (taking 2% doping), and has 6 arms 
extending along [100], [010] and [001]. The arms along the 
z-axis are slightly longer than those along x and y, but 
on the scale of Fig. [3] this is not a perceptible difference. 
In contrast, Mattheiss's Fermi-surfaces (Fig. 6 of Ref. l40|) 
have 4 arms along [100] and [010] and none along [001]. 
Gregory et a?P^ studied samples with electron density 
6 • 10 18 cm" 3 , corresponding to x = 3.6 • 10~ 4 . Due to the 
large crystal field splitting, the Fermi level in Mattheiss's 
calculation is then still below the second band. Yet Gre- 
gory et al. observed low-frequency quantum oscillations 
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with frequencies 40 Tesla and 45 Tesla. The weak field- 
orientation dependence indicated that these are associ- 
ated with rather isotropic Fermi surfaces, which they 
associated with the light-electron band. To have this 
band occupied they postulated that Mattheiss' estimate 
of the splitting of the two lowest bands introduced by 
the tetragonal distortion needed to be revised downward. 
Looking now at our calculation we notice that, since 
x = 3.6 ■ 10~ 4 > x c i, the light-electron band is indeed 
occupied for this doping level. As shown in Fig. [5] the 
diameter of the second Fermi surface is practically in- 
dependent of direction for this low doping range with 
a radius k — 0.134/a = 3.54 • 10 6 cm' 1 and extremal 
area A — 3.95 • 10 13 cm' 2 . Using the Onsager relation 
F = Ah/(2ne) the corresponding quantum oscillation 
frequency is 41 Tesla. Since these samples consisted of 
many domains with different orientation of the tetrag- 
onal axis, a doublet due to the anisotropy is expected 
and observed. It thus appears, that the new ab initio 
band structure settles an old conundrum regarding the 
quantum oscillations of n-type SrTi03. 



ujp.bj refers to the fact that, since the LDA-calculation 
does not take into-account electron-phonon interaction, 




FIG. 5: Enlarged view of Fig. [3] indicating the position of 
the Fermi level for x = 3.6 ■ 10~ charge carriers. The corre- 
sponding value ka = 0.134 is in excellent agreement with the 
hitherto unexplained de Haas-van Alphen frequency reported 
by Gregory et aZp^. 



IV. MASS RENORMALIZATION INDUCED BY 
ELECTRON-PHONON COUPLING 



In Ref. [M] we compared the Drude spectral weight to 
the same quantity calculated using LDA. The expression 
for the spectral weight along the Xj-axis is 
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where the sum is over momentum, band-index and spin 
and /(e) is the Fermi-Dirac distribution. The index b in 



it calculates the bare mass. The ratio . 
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id 2 e j is the experimental Drude spectral weight, then 
corresponds to the mass renormalization factor m* /rnb- 
This procedure was followed in Ref. [3H Since the lin- 
ear term of the specific heat is a direct measure of the 
density of states at the Fermi energy, 7 = B £ N F , the 
ratio 7e/7fc of the experimental over the LDA value pro- 
vides a second way to measure the mass enhancement. In 
Fig. [6] the LDA-calculat ion of the DOS at e F is plotted 
as a function of doping, together with values obtained 
from experimental specific heat data. Clearly the DOS 
as given by experiments is about a factor 2 higher than 
the LDA-prediction. The corresponding mass enhance- 
ment together with the results of the other two methods 
are summarized in Fig. [6] The verdict is clear: There is 
a factor of 2 to 3 mass enhancement with a tendency to 
become smaller for higher doping. Electron-phonon cou- 
pling is the only plausible suspect for the enhancement. 
Indeed, recent calculations confirm this^: Based on the 
Frohlich interaction the essential characteristics of the 
observed optical conductivity spectra of SrTii-ajNb^Oa, 
in particular intensity, lineshape and energy of a peak at 
130 meV, was explained without any adjustment of mate- 
rial parameters. The electron-phonon coupling coupling 
constant was found to be of intermediate strength. 



For the correct understanding of the peculiar temper- 
ature dependence it is important to find out whether 
or not the charge carriers are to a good approximation 
described by Bloch waves. This corresponds to the re- 
quirement that the mean-free path at the Fermi-level, 
I = v* f t is much bigger than the Fermi-wavelenth, in 
other words v* f t >• 2n/k F . Multiplying both sides of 
the expression with k F /2 we obtain e F h~ 1 T 3> tt, where 



?/v F — nib/m*. In the previous section we 



have obtained the doping dependence of e F . Combining 
this with the m* of Fig. [6] and h/r of Fig. [2] we obtain k F l 
as a function of temperature for different dopings, shown 
in Fig. [7] We see from this graph that at 4 K the elec- 
trons are strongly Bloch-like. At low temperatures the 
largest k F l occurs for the lowest carrier concentration. 
This is expected since in these samples the number of 
charge carriers is equal to the number of Nb-ions, which 
act both as donor atoms and scattering potentials. The 
opposite trend occurs above the isosbectic point (24 K, 
k F l — 33). While at 4 K we obtain high values of k F l in 
the range from 40 to 150, at 100 K we have k F l of the 
order 2ir implying localization in Fermi-wavelength sized 
wavepackets. In the low temperature range it is there- 
fore reasonable to extract interaction parameters from 
the coefficients of the T 2 dependence of 1/r. Above ap- 
proximately 100 K the material enters into a regime of 
incoherent transport. We therefore restrict the analysis 
of 1/r in section [V] to temperatures below 100 K. 
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FIG. 6: Top panel: Doping dependence of the density of states 
(DOS) at the Fermi energy. The solid curve corresponds to 
the tight-binding bandstructure fitted to the Wien2k ab initio 
bandstructure, with p arameters of the first row of table |IV| 
Squares^, circles^^, pentagon^ Density of states obtained 
from the linear term in the specific heat. Middle panel: Dop- 
ing dependence of the Drude spectral weight, Up. The solid 
curve corresponds to the bandstructure results. Diamonds 
are the experimental values^. Bottom panel: (i) Ratio of ex- 
perimental DOS over bandstructure DOS (experimental and 
theoretical values taken from top panel, the meaning of the 
symbols is the same), (ii) Ratio of bandstructure over ex- 
perimental ujp (values taken from middle panel), (iii) Ratio 
of bare and experimental (dressed) Fermi-velocity, VF,b/vF,e 
(triangle, data from Ref. 37). The grey curve is a smooth 
interpolation, m* /m& = 2.0 + 1.2 exp (— a/0.005). 



T RELAXATION RATE AND TWO-BODY 



INTERACTIONS 



The low carrier density leads to a situation where the 
kinetic energy of the charge carriers is slaved to the 
relevant vibrational energy scale. The usual Migdal- 
Eliashberg expansion in the electron-phonon coupling 
constant is therefore not applicable. A different approach 
is required whereby in first instance the electron-phonon 




FIG. 7: The quantity kpl where I is the mean free path, 
calculated using k F l = 2e* F h~ 1 T, and using the experimental 
r of Fig. [2] and the calculated Fermi energies corrected by the 
mass renormalization of Fig. [6] The high values of k F l below 
100 K imply the itinerant character of the charge carriers. 



coupling is treated for each individual electron, resulting 
in charge carriers renormalized by electron-phonon cou- 
pling, which condense into a Fermi-liquid of "polarons". 
The polarons interact with each other via the Coulomb 
interaction and by virtual exchange of phonons. With re- 
gards to optical phonons the hierarchy of energy scales is 
inverted as compared to the situation in common metals, 
in that Hoj > e* F where u) is the optical phonon frequency 
and e* F the Fermi energy of the polarons. While it is clear 
that the Migdal-Eliashberg expansion can not be used, 
the solution of the many-body problem in this limit is a 
complicated problem which we will not attempt to solve 
here. Instead we turn the problem around and anticipate 
that the correct solution should share certain properties 
in common with the problem of interacting composite 
fermions such as 3 He. The essential properties should 
then be those of fermions interacting through some effec- 
tive interaction mediated by the optical phonons, which 
on the scale of e* F can be considered effectively non- 
retarded. An immediate consequence is then a T 2 con- 
tribution to the inelastic relaxation rate, resulti ng from 

271! 



the phonon-mediated fermion-fermion interactio n l 28 ^ 
The situation in SrTii-^Nba-Os is more complicated 
than in 3 He in that the Fermi surface is crossed by 3 
bands for x > 2.6 • 10~ 3 and by 2 bands for 4.0 • 10~ 5 < 
x < 2.6 • 10~ 3 . However, we assume that the trans- 
port and superconducting behavior are dominated by the 
most highly occupied band, and use the single-band ex- 
pressions in Appendix [B] to derive the effective coupling 
constants. 

If indeed the temperature dependence of the relaxation 
rate is a manifestation of fermion-fermion scattering, it 
should be possible to obtain from it a parameter char- 
acterizing the interaction strength. The formalism has 
been elaborated in the context of the observation of van 
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Kempen et al. of a T 2 contribution in the resistivity of 
potassiurrPS. In particular, the expression for the relax- 
ation rate is (see appendix [B]) 



h 

T 

0,2 



a 2 T 2 



TT 3 k 2 



= A 



(3) 



The parameter u < 1 describes the fraction of the mo- 
mentum changes which is transferred to the ionic lattice. 

The dimensionless parameter A T represents the inter- 
action effective in polaron-polaron scattering. Since we 
have determined a 2 and e* F in the previous sections, we 
are ready to calculate A 2 u using this expression. The 
result is shown in the lower panel of Fig. [8] The trend 
of At-u 1 / 2 going to zero in the zero-doping limit may be 
a consequence of either A r or u diminishing in the low 
doping limit, or a combination of these two. As pointed 
out in AppcndixjB] several factors make u ^ 0: (i) Baber 
scattering involving heavy and light electrons (ii) Umk- 
lapp scattering, and/or (iii) disorder scattering by donor 
atoms. Since for u — > (i) the mass-anisotropy (section 
|III[ ), (ii) the probability of intra-pocket Umklapp, and 
(iii) the impurity scattering 1/tq (see Table Ref. [T]) all 
vanish, we may expect u — > in this limit. Since we 
will see in the following section that the polaron-polaron 
interaction as calculated from T c is almost independent 
of doping, we tentatively attribute the observed doping 
dependence of A 2 w to a suppression of u for low doping. 



VI. SUPERCONDUCTIVITY 

Superconductivity is observed in n-type SrTiC>3, with 
a dome shaped T c between and 0.02 electrons per 
SrTiC>3 formula unifr^, with maximum values of about 
0.7 K when doped with NrM This doping dependence 
and the T c itself are relatively robust features of the 
doped 3-dimensional bulk materials as well as the 2- 
dimensional SrTiC>3 /LaA103 interfaced^. Superconduc- 
tivity in doped bulk SrTiOs has been anticipated by 
M. L. Cohen on the basis of an attractive electron- 
electron interaction arising from the exchange of intraval- 
ley and intervalley phononiP^ and motivated by early 
bandcalculations indicating a many-valley bandstruc- 
ture in SrTiOs. The intervalley mechanism has been 
further elaborated in the context of SrTiOs in a num- 
ber of papers^SHlSl However, over the years evidence 
has been accumulating that all bands are at the cen- 
ter of the Brillouin zone. Several alternative mechanism 
not involving a multi- valley bandstructure have been pro- 
posed, to mention a few: (i) J. Appel 51 noticed that the 
Brillouin-zone folding associated with the tetragonal dis- 
tortion creates two bands of zone-folded optical phonons 
with a quasi-accoustic dispersion at the zone center. One 
of these bands has a finite matrix element for intraval- 
ley scattering and can therefor in principle mediate su- 



perconducting pairing, (ii) Z. Zinamon, while maintain- 
ing Appel's idea of soft phonon exchange, argued that 
the relevant charge carriers are small polarons, and pro- 
posed a theoretical model relevant for this limilP^l, (iii) 
T. Jarlborg has demonstrated by electronic structure cal- 
culations that the electron-phonon coupling is enhanced 
for long-wavelength phonon mode despite the low den- 
sity of states, which is consistent with the appearance 
of superconductivity at low doping^. We see, that the 
mechanism for pairing in this material is far from clear. 
We therefor adopt here a phenomenological approach 
whereby we deduct the coupling constant characteriz- 
ing pairing interaction from the experimental T c 's and 
compare it to the coupling constant obtained from the 
transport data. While this approach does not solve the 
question as to the exact nature of the phonon-mediated 
interaction, it does allow to establish whether supercon- 
ductivity and transport properties can be treated in a 
unified approach of an interacting Fermi-liquid. 

One of the consequences of the Fermi energy being 
smaller than the relevant phonon energy uj c , is that the 
energy cutoff of the pairing interaction is given by e* F on 
the occupied side of the Fermi level, while it is given by 
lo c on the unoccupied side. This introduces a dependence 
of T c on e F which may in part be responsible for the de- 
crease of T c for x — > 0. A more detailed discussion of the 
consequences in the limit of weak coupling (T c <C Tp, as 
is the case in all samples that we discuss here) is provided 
in Appendix [C] Before we set out to discuss this, it is 
important to establish whether the superconductivity is 
closer to the BCS limit, or to the limit of Bose-Einstein 
condensation of bipolarons. The latter has been pro- 
posed for low carrier concentrations (n < 10 18 cm ) in 
Zr doped SrTiOgEH and for the high Tc cupratesP. One 
way to investigate this question is, by estimating with 
how many other pairs each Cooper-pair overlaps. In the 
BEC-limit there is essentially no overlap, whereas in a 
BCS superconductor it is given by the volume occupied 
by a pair divided by the available volume. The former is 
just (47r/3)£o where £o = Iivf/tt^o, and the latter is 2/n 
where n is the electron density. Using standard relations 
between density and Fermi energy we obtain 



V(occupied) 4 
V (available) ~ 9tt 4 \A 



(4) 



Binnig et oZP have observed gap values close to 
Ao/ksTc = 1.76 in their tunneling spectra for exactly 
the data in Fig. [HJ hence we can use this substitution for 
A . The result is shown in the middle panel of Fig. [8] 
in a broad doping range using data collected by Koonce 
et aZP and Binnig et alW. Ipso facto each Cooper-pair 
overlaps with 10 5 to 10 10 others, which places these su- 
perconductors clearly outside the realm of Bose-Einstein 
condensation for the range of carrier concentrations con- 



sidered here. Substituting in Eq. C5 of Appendix O the 



values of T r and the value of e* F discussed in section IV 



we calculate the corresponding coupling constant A for 
the pairing interaction. Since we lack certainty about the 



n (cm" ) 



lxlO 20 2xl0 2 ° 



3xl0 2 ° 4xl0 2 ° 



0.6 

& 0.4 
H 

0.2 
0.0 

io 10 

c a i0 8 
1 10 7 

> io 6 
io 5 

0.4 
g 0.3 
0.2 
0.1 




Q 



i i i i | i i i i | i i i i | i i i i | i i i i 



I I I I I I I I I I I I I I I 



I I I I I I I 




0.4 
0.3 
0.2 
0.1 
0.0 



0.0 

0.000 0.005 0.010 0.015 0.020 0.025 
x 

FIG. 8: Top panel: T c of a large number of samples. Pen- 
tagons: data by Binnig et a/P. Circles: data by Koonce et 
Middle panel: Coherence volume times the density of 
Cooper pairs. The high value indicates that each pair is over- 
lapping with a huge number of other pairs. Lower panel: Dop- 
ing dependence of the coupling constants X T obtained from 
the T 2 term of the resistivity and Ao from T c . Circles: A T 
using a,2 listed in the second column of table [XT] and Eq. [3] Di- 
amonds: idem using a 2 . Bars: coupling constants A from T c 
at the corresponding carrier concentration of SrTii-^Nb^Os 
using the T c 's of the top panel and Eq. |C5| e* F = eFmb/m* 
using tF from Fig. [3] and m* /mb from Fig. |6| Upper (lower) 
limits of the bars indicate the value for g obtained with lo c = 
5 meV (80 meV). 



nature and frequency of the phonons causing the pairing 
interaction, we have substituted two extremal values for 
the vibrational cutoff-energy in the gap equation: co c =5 
meV and 80 meV. The resulting uncertainty of Ao is not 
very large; for all dopings we find 0.1 < A < 0.2 with 
negligible doping dependence. The results are shown in 
Fig. [sj together with XrU 1 ^ 2 . The different doping depen- 
dence of Ao and A^-u 1 / 2 has a simple explanation in that 
we expect the parameter u to vanish for x — > 0. 



VII. LANDAU PARAMETERS 

In principle we want to determine the full set of rel- 
evant Landau parameters, either in the form A\ or as 
Ff . However, even while we have assumed that the only 
relevant angular momentum values are I = 0,1, there 
are still 4 parameters while until now we have deter- 
mined two quantities which depend on them, namely Ao 
and A,-?/ 1 / 2 . The first question concerns the symmetry 
of the pairing itself: The expression relating the super- 
conducting coupling constant to the Landau parameters 
are different for singlet and triplet pairing, so one has to 
make a choice as to whether one assumes triplet or singlet 
superconductivity. Triplet pairing can be excluded be- 
cause the only available mechanism in the present case is 
electron-phonon coupling. The second question concerns 
the value of u. In the previous section we attributed 
the suppression of Ai-u 1 / 2 as x — > to the suppression 
of momentum transfer to the ionic lattice. Vice versa, 
in Fig. [8] we see that \ T v}l 2 has saturated for x > 0.2. 
Accuracy by which the Umklapp fraction can be calcu- 
lated is probably within a factor of 2, even in the alkali 
metals which are relati vely simple due to the nearly free 
electron characte r^ 2 * 56 !. We make the simplest possible 
assumption that u»l, implying that for the higher dop- 
ing levels A T ~ 0.4. We will base the analysis of the 
Landau parameters on this value. With 4 parameters to 
determine and 2 experimental constraints we need two 
additional pieces of information. One of them is sup- 
plied by the sum rule for the Landau parameters derived 
by Brinkman, Platzman and Rice 5 - (BPR sum rule) for 
charged fermions, which for the sp-model implies 



Al+A% + A1 = -1 



(5) 



We are still one constraint short. One might hope to find 
such a constraint in, for example, the mass-enhancement 
measured with specific heat. The problem is however 
that one needs to compare the values of the electronic 
specific heat with and without polaron-polaron interac- 
tions. Since the mass of a polaron is already enhanced 
compared to the bare band mass by a factor of 2 approx- 
imately, to extract the contribution of polaron-polaron 
interactions, especially if it is much smaller than 1 as 
it turns out to be the case here, is difficult and at the 
present state of affairs not feasible. We therefore calcu- 
lated Aq, Af, and A1 as a function of Aq while fixing 
the constraints imposed by Ao = 0.15 through Eq. C4 
by A r = 0.4 through Eq. |B7| and the sum rule Eq]5 
The parameter Aq is varied in the range of positive mass 
enhancement (ra*/m — 1 = F{/3 > 0), and positive com- 
pressibility (k = K},(1 — Ag)m* jm > 0, where Kf, is the 
bare value). Since we assume that the pairing symmetry 
is of the singlet variety, we only consider the solutions 
for Ao > Ai. The result shown in Fig. [9] allows to deter- 
mine all parameters once the value of Aq has been set. 
The difference A\ — Aq represents an exchange interac- 
tion which tends to align spins parallel for positive values. 
Its value increases for Aq —> —1.27. Correspondingly, for 
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— 1.27 < Aq < —1.15 an alternative set of solutions is 
obtained corresponding to triplet pairing (not displayed 
in the figure) which, as already pointed out above, we 
reject on theoretical grounds. 

Even with this broad range of possibilities for Aq al- 
lowed by the experimental constraints, the windows for 
Aq, A\, and A\ are limited. In table III we compare all 
parameters discussed here to the case of 3 He at ambi- 
ent pressure. We see, that the values for SrTiC>3 are of 
the same order, but smaller than in liquid 3 He. The fact 
that we obtain "reasonable", i.e. not excessively large 
or small numbers, of the Landau parameters, gives fur- 
ther support to the notion that that n-type SrTi03 is a 
Landau Fermi liquid, and superconductivity and the T 2 
resistivity in this compound have a common origin. 



o 



- 




Parameter 


SrTi0 3 


3 He,l atm. 




singlet 


triplet 


AS 


{-1.27; 1.0} 


0.91 


Af 


0.45 ±0.25 


2.0 


AS 


-0.67 ±0.22 


-2.03 


Af 


-0.62 ±0.28 


-0.55 



TABLE III: Ranges of the Landau parameters of SrTiOs al- 
lowed by the constraints imposed by the experimental data 
combined with the BPR sum rule (2d column, this work) and 
values of the same of He (3d column, see Ref. 1581 original 
data in Refs. I59I60|) . 



VIII. CONCLUSIONS 

We have performed a detailed analysis of the T 2 behav- 
ior of the resistivity of n-type SrTiC>3. Novel bandstruc- 
ture data are presented, and it is shown that the new 
band structure solves an old conundrum of the de Haas- 
van Alphen frequencies dating from 1977. The mass, den- 
sity of states, and plasma frequency of the quasiparticlcs 
are found to be renormalized by approximately a factor 
2 due to electron-phonon interaction. The quasiparticles 
turn out to be in the anti-adiabatic limit with respect to 
electron-phonon interaction with a quasi-instanteneous 
interaction mediated through phonons. Analysis of the 
T resistivity and T c provides values of the Landau pa- 
rameters of n-type SrTiC>3 which are comparable in size 
to those of liquid 3 He. 
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FIG. 9: Solutions for Af, Aq, Af, the triplet superconducting 
coupling constant and the mass-enhancement as a function 
of Aq when taking the experimentally determined A T = 0.38, 
Ao = 0.15 and the BPR sum rule as constraints assuming 
singlet pairing for the ground state. 



Appendix A: Bandstructure 

Calculations have been performed in the high- 
temperature perovskite (HTP) structure (a =3.905 A) as 
well as in the low-temperature tetragonal (LTT) struc- 
ture (group no. 140, I4/mcm, a =5.529 A, c = 7.824 
A, O x — 0.244); optimizing the O position in the cal- 
culations yields O x = 0.223, indicating, not surprisingly, 
that even zero-point fluctuations substantially reduce the 
average distortion. Nb doping was simulated in the vir- 
tual crystal approximation, changing the nuclear charge 
of Ti from 22 to 22 + x, or that of Sr from 38 to 38 + x 
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(the results did not change, proving that this is a good 
approximation in the considered range of dopings). 

Most calculations were performed with RK maK — 7 for 
the wave-function expansion, and i?G max = 14 for the 
charge density expansion; calculations with i?i"T max = 8 
and i?G max = 16 did not show a discernable difference. 
Fermi-surface integrals were evaluated using the k-point 
meshes up to 28 x 28 x 28. The plasma frequencies were 
evaluated as the Fermi-surface averages of the squared 
Fermi velocities. The Fermi velocities were calculated 
using the WIEN2k optics package; accuracy of the Fermi 
velocities was tested (in the cubic case) by numerical dif- 
ferentiation of the energy eigenvalues. 

For many calculations of physical properties it is useful 
to have an integration in momentum space which is rapid. 
We therefore used a parametrization of the ab initio band 
structure of the ti g manifold in the region around the 
zone center, based on the following tight-binding model 
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sin 2 








k.l 











£ fe,2 














ha 
~2~ 



U s sin 



! (ZD £ £ 
with the corresponding parameters given in table |IV| 



(Al) 



Source 


ts 




f 


D 




meV 


meV 


meV 


meV 


Wien2k code 


35 


615 


18.8 


2.2 


Mattheiss 


35 


500 


28 


-33 



TABLE IV: Tight-binding parameters describing the disper- 
sion of the t2 g bands near the P-point of SrTiC>3. First 
row: Parameters fitted to the Wien2k code bandstructure 
presented here. Second row: Parameters fitted to Matheiss' 
calculations^ for the tetragonal phase with a tilt angle of 
2.1°. 



Appendix B: Landau parameters 

In the Landau-Fermi liquid theory of interacting 
fermions, the bare interaction is expressed in terms of 
the dimensionless Landau parameters Fj 1 where the in- 
dex / indicates the angular momentum and j the parity 
of the scattering process. The mass renormalization de- 
pends only on 



- 1 = 



Ft 



(Bl) 



The specific heat is renormalized by the same factor 
as the effective mass. The interaction between dressed 
quasiparticles 



Ff 



l + if/(2Z + l) 



(B2) 



is important for scattering between quasiparticles and 
constitutes the pairing interaction for superconductivity. 
Following Dy and Pethic#3 we limit the scattering pro- 
cessus to the I — and I = 1 values in the spherical ex- 
pansion of the scattering amplitudes (the so-called s — p 
approximation) 



M 



i 

1 

Np 



(A s - 3Ag) + (Al - 3A1) cos6] (B3) 



(Al + A") cos 9] coscj> 



where N F — m*k F /(Tr 2 H 2 ) is the density of states at the 
Fermi level and the angles 9 and <f> represent the kinemat- 
ics of the scattering events. The relevant quantity in the 
theory of ine lastic scattering is the transition probability 



W(9^) = -[A t 



(9,<f>) + A t (9,<j>)Y 



2h 



-At 



(B4) 



Due to collisions between quasiparticles, the relaxation 
rate of the dressed quasiparticles has a T 2 temperature 
dependence, which is the most characteristic property of 
a Landau Fermi liquid. We will f ollow here the approach 
of Lawrence and Wilkins ^ 7 * 28 * 45 !. These authors define 
a surface-averaged relaxation rate for an electron at the 
Fermi surface^, and derive the relation between r and 
W(9,<t>) 



1 = (m* f(k B T) 2 u / W(0,<l>) \ 
r 12ir 2 h e \ cos (0/2)/ 



(B5) 



where the dimensionless coefficient u < 1 represents the 
efficiency of momentum transfer to the ionic lattice of 
the relaxation process. In a translationally invariant sys- 
tem of interacting electrons u = 0, because the current 
operator commutes with the Hamiltonian of such a sys- 
tem. However, the fact that a solid does not possess 
full translation symmetry has important consequences. 
Already in 1937 Baber demonstrated a mechanism for 
finite resistivity in a two-band model in which s elec- 
trons are scattered from heavier d holes by a screened 
Coulomb interaction^. The Baber mechanism works 
more generally for a system of light and heavy electrons 
with the heavy particles acting as momentum sinks, and, 
as pointed out Giamarchi and Shastry, similar results 
are expected for typical, noncircular bands^. This last 
point is relevant to the case of SrTiC>3 in view of the 
strong fc-dependence of the mass across the Fermi-surface 
(see Fig. Eh. In single band Umklapp processes allow 
momentum transfer to the crystal coordinate system^]. 
Likewise, the potential landscape caused by impurities 
(for example the donor and/or acceptor atoms in doped 
semiconductors) provides a channel by which momentum 
gets transferred to the ionic lattice in electron-electron 
collisions. Refs. 27 28 45] concentrated on alkali-metals, 
and therefor did not take into account Baber scattering. 
They used the symbol A for the Umklapp fraction. To 
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avoid confusion with the superconducting gap, we indi- 
cate here the fraction of momentum in electron-electron 
collisions transferred to the ionic lattice due to Umklapp 
and other mechanisms with the character u. Integra- 



tion of the angular integrals in Eq. B5 is straightforward 
though tedious, with the result 



w(0,4>) \ 



where 



I2\i 



cos 



24 



(0/2)/ 



12A 



>h 5 



(At) 



7 



(B6) 



8 



12 



7 

A s A s — - A a A a 



4Q 

f 2 ±(a«Y-\a* a* 

1 AS AO, I _ AO, AS 



(B7) 



If A\ is the only non-zero parameter we obtain A 7 /(l + 
A 7 ) = 4/3 A = 4 V /277A T . Finally, by substituting 
Eq. |B6| in Eq. |B5| we obtain 



h 

T 



-\lu{irk B Tf 



(B8) 



This is the central expression enabling extraction of \ 2 u 
from the experimental values of the amplitude of the T 2 
term in the resistivity. 



Appendix C: T c equation 

The gap equation for an isotropic gap is 

'(e-rf + A 2 



r d£ _ N(e)V(e) 



"2V(e-/i) 2 + 



A 2 



: tanh 



2k B T 



(CI) 



The chemical potential /i(T c ) is to be determined at the 
critical temperature by adjusting it such as to fix the 
number of electrons 



N(e) 



1 



1 _|_ g/3 c (e-M) 



de 



N(e)de 



As it turns out to be the case for the data considered in 
the present paper, ksT c << e* F ; consequently the output 
of the self-consistent solution is [i(T c ) w e^. The critical 
temperature is obtained by solving the gap equation for 
A = 0. In the present case the bottom of the band 



constitutes the lower limit of the integral over the density 
of states. We define it as the zero of energy, so that 



deN(e)V(e 



tanh(fe(6-^)/2) 
2(e - (,) 



de 



(C2) 



where fcsTc = 1//3 C . The usual approximation for the 
retarded interaction consists of substituting N(e)V(e) = 
A for | e — /j. | < u) c where uj c is cutoff energy of the pairing 
interaction, and taking A = for |e — fi\ > uj c . The 
expression for T c is then 



/cflTcj- = 1.13cj c exp 



-1 



(C3) 



where uj c is the cutoff energy of the pairing interaction, 
and the coupling constants for the I = (singlet) and 
I = 1 (triplet) pairing channels are^ 

Ao = \ {A\-Al) + \(A«-A1) 

Ai = ^{At-Al)-±{A a -AD (C4) 

In the case of n-type SrTi03 an interesting asymmetry 
is introduced by the condition that the energy scale of 
the phonons mediating the interaction is in the anti- 
adiabatic limit: we have fj, < uj c . Since on the occupied 
side N(e) = for — lj c < e — fi < —fi, the region of finite 
N(e)V(e) is limited to — fi < e — /i < lj c . Another aspect 
to take into account is that N(e) — c^/e. We therefore de- 
fine the dimensionless coupling constant A at jj, through 
the relation N(e)V(e) — \\[Tf\± and we make a transfor- 
mation of variables x = f3 c (e — (i). The equation for T c 
then becomes 



1 

A 



1 



0eUJ ° I ^tanh(x/2) J 



2.r 



(C5) 



which has the following solution in the weak coupling 
limit (A < 1) 



T c = 0.612^ exp 



exp 



-1 



(C6) 



The values for Ao shown in Fig. [H] were obtained by solv- 
ing Eq. C5 numerically, and agree within 3% accuracy 
with the weak coupling expression Eq. |C6| for oj c =80 
meV. 
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